Optimizing the Multi-Photon Absorption Properties of NOON States 
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In this paper we examine the A-photon absorption properties of "NOON" states, a subclass of 
path entangled number states. We consider two cases. The first involves the A*'-photon absorption 
properties of the ideal NOON state, one that does not include spectral information. We study how 
the A'^-photon absorption probability of this state scales with A'^. We compare this to the absorption 
probability of various other states. The second case is that of two-photon absorption for an A = 2 
NOON state generated from a type II spontaneous down conversion event. In this situation we 
find that the absorption probability is both better than analogous coherent light (due to frequency 
entanglement) and highly dependent on the optical setup. We show that the poor production 
rates of quantum states of light may be partially mitigated by adjusting the spectral parameters 
to improve their two-photon absorption rates. This work has application to quantum imaging, 
particularly quantum lithography, where the N-photon absorbing process in the lithographic resist 
must be optimized for practical applications. 

PACS numbers: 42.50.St, 42.50.Hz, 85.40.Hp, 42.50.Dv 
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I. INTRODUCTION 

We investigate the multiphoton absorption probabili- 
ties of maximally path entangled number states - also 
called NOON states, after the way the state vector is 
written: |7V :: 0) = {\N,0) + \0,N))/V2. The two posi- 
tions in each state vector represent two spatial modes in 
an optical interferometer. So, for example, for modes 1 
and 2 we abbreviate |a;)i|y)2 as \x,y). These states are 
of interest due to the fact that they greatly improve the 
resolution and sensitivity of interferometry for metrology. 
Also, it has been shown that they would improve the res- 
olution with which lithographic features may be written. 
We show that (monochromatic) NOON state absorption 
fares poorly as N increases. Thus when considering pos- 
sible applications of these states, such as to quantum 
lithography P, 0] or metrology [1, Q , we need to keep an 
eye towards maximizing these absorption rates by vary- 
ing their spectral parameters. Wc do this knowing that 
it has been found that squeezed light can exhibit novel 
two photon absorption properties, such as linear growth 
of absorption rate with intensity and decreasing absorp- 
tion, for increasing field [5j, iS, 13] ■ 

We consider in detail the case of a |2 :: 0) state used 
in a quantum lithography or quantum metrology setup. 
We include spectral information and derive a general ex- 
pression for the two-photon absorption probability. Then 
we numerically maximize the probability function and 
find the setup which maximizes the two-photon absorp- 
tion. The absorption probability can be improved by 
several orders of magnitude by carefully adjusting the 
filter bandwidths, pump pulse length, and the length of 
the crystal. We go beyond most previous studies in that 
we obtain the two photon absorption probability directly, 
instead of only considering the second-order correlation 
function. 

In Section II we compare the absorption properties of 



ideal NOON states to other states of light. In section 
III wc calculate the biphoton amplitude in the general 
case and then examine the absorption properties of this 
type of light in two regimes: the pulse-pumped, and the 
continuous-wave-pumped. Some of the more lengthy cal- 
culations are left to the appendices. 



II. ABSORPTION PROPERTIES OF IDEAL 
NOON STATES 

Initially one may consider the ideal NOON state: :: 
0) [l[. This state contains no spectral information. It 
is an abstraction which can only exist in an optical cav- 
ity. Nonetheless it will provide some insight into how the 
absorption properties of NOON states compare to other 
sources. 

Agarwal studied how multiphoton absorption rates are 
influenced by the specific properties of the incident light 
0. He found that the equation of motion for the field 
can be written as 



dt 



-2riA(")(a^"a"). 



Here, A*-"^ is the absorption coefficient for the n-photon 
absorption process and contains information about the 
medium that is acting as the absorber. The rate of 
change of the number of photons in the field is propor- 
tional to the absorption rate. So the probability of an 
71-photon absorption event occurring is, 



(at a)" ' 



where k is some constant which we set to one in the inter- 
est of simplicity. This represents roughly the probability 
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FIG. 1: Multiphoton absorption probability of various sources 
as a function of the average number of photons in the field 
(A'"*'' order absorption vs. N average photons). All probabil- 
ities normalized relative to coherent states (solid line). The 
dashed line is a thermal state. The dotted line is a Fock 
(number) state, and the dot-dashed line is a NOON state. 

of n photons in a light field being absorbed by an n- 
photon resist. We can use this information to produce a 
graph to see how the multi-photon absorption probability 
of NOON states scale with N, when compared to quan- 
tum states that are not path entangled such as thermal, 
number, and coherent. 

Thermal states are described by the following density 
matrix 



/5the™al = ^^e-^^/'=^|j)(j|, 
g~hLj/2kT 2^ 

Obtaining the matrices for number, coherent, and NOON 
states is straightforward. For a two mode state the anni- 
hilation operator is given as a = ai -f 02 where one and 
two label the two paths the photon may take 

See Fig. [1] Thermal states clearly have the greatest 
rates of multiphoton absorption. This can be attributed 
to the fact that thermal states exhibit bunching. That is 
that photons from thermal radiation tend to be tightly 
correlated in time. Fock (number) states fare the worst. 
This feature of Fock states is connected to the fact that 
number states represent standing waves where the loca- 
tions of the individual photons are evenly spaced out 
(or anti-bunched) in space and time with no definable 
phase. The multiphoton absorption properties of coher- 
ent states stay constant with respect to photon number. 
Since photons in coherent states are randomly dispersed 
in space and time the chance of two or more photons 
being correlated is simply proportional to the intensity 
(average photon number). Since we are considering how 
ri-photon absorption scales against average photon num- 



ber, the graph is flat, providing a convenient measuring 
stick to gauge other fields. 

For NOON states the multiphoton absorption behaves 

as 

Pjv = 1 for = 1, 

A^! 

= 2^rr for N >2, 

where k has been set to one. NOON states fare a fac- 
tor of two better than Fock states, although absorption 
rates are still far from optimal. The reason NOON sates 
have this factor of two is due to the path entanglement. 
Mathematically this two comes from the normalization 
constant that path entanglement requires. 

These results seem to reinforce the findings of Tsang 
in Ref . [13 : generally the absorption properties of NOON 
states are poor. It should be re-emphasized that these 
states, containing no spectral or temporal information, 
are idealizations. Fig. [T]can only be seen as providing a 
rough idea of how absorption scales. 

Quantum lithography or metrology will only be useful 
if the detector or material needs to be exposed to the 
field for a reasonable period of time. The above ideal- 
field results seem to make this unlikely. However, they 
show only that the quantum mechanical properties (i.e. 
the bare state vector) of NOON states lead to poor ab- 
sorption rates. They say nothing about how the spectral 
properties of realistic NOON state pulses effect the multi- 
photon absorption probability. We thus examine in de- 
tail a specific well known case: the |2 :: 0) state. Though 
this state is not practical for metrology or lithography 
its optimization would provide a proof of principle that 
the absorption properties of higher N states could be 
improved in a similar manner. 

III. ABSORPTION PROPERTIES OF 
REALISTIC |2 :: 0) STATES 

It is possible to write realistic states, which include the 
spectral information of the light of interest. Furthermore, 
after these states are used to obtain absorption rates, the 
arrangement of optical elements which optimizes absorp- 
tion can be found. 

Several works have examined how biphotons pro- 
duced by parametric down conversion or electromagnet- 
ically induced transparency may be compressed or oth- 
erwise modified so that they exhibit tighter correlations 
[O, E, [H, [13, E [ii- There is an excellent paper by 
Dayan which studies the properties of a semi-stationary, 
undepleted beam of squeezed light produced by a spec- 
trally narrow pump . 

We consider the case of a |2 :: 0) state used in a 
quantum lithography [l], Q or remote quantum metrology 
[i,[3l setup. The |2 :: 0) state we investigate is produced 
by co-linear type II degenerate down conversion and a 
beam splitter (sec Fig. [5]). This setup is very simple, but 
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FIG. 2: The basic setup. A nonlinear crystal (BBO in this 
case) creates a degenerate pair of photons. Each photon is 
subjected to a filter. A polarization rotator ensures that the 
two photons are indistinguishable. A beam splitter creates a 
1 2 :: 0) state that results in an interference pattern. 



by tuning these few basic optical elements we can see a 
large improvement in two photon absorption rates, with- 
out recourse to exotic techniques. We make no assump- 
tions about stationarity or about the relative sizes of the 
field bandwidths (apart from one very broadly applicable 
assumption - that the field's bandwidth is narrower than 
the atomic transition frequencies; this will be discussed 
fm'thcr in Appendix B). 



A. The Type-II Biphoton 

The output state of the crystal during type-ff down 
conversion is described by [3 



J J duj'dujg*{uj') 
X S*'^-* {uj f ~ uj\ uj' \ uj f ~ oj , uj)g{uj). 



(2) 



Where g{uj) is the atomic response function, tof is the 
frequency of the final state, and 5''^^) is the spectral cor- 
relation function (the fourier transform of the temporal 
correlation function), which in our case is 



(3) 



where, 



Here, lu[ and represent the negative frequency compo- 
nents of the field associated with t['^ and t'2'^. The factors 
LUi and UJ2 represent the positive frequency components of 
the field associated with ti'^ and t2'^. The above equation 
does not include the effects of natural linewidth, which 
will be discussed later. We should note that other equa- 
tions, also derived by MoUow [23|- assume that the field is 
stationary, something which is not true in general for our 
calculation here. A is the biphoton amplitude, defined 
below. 

We start with the two-photon correlation function [2^ 



|^')=CV / dujp f dze ^(^^^) e''^""' 
kk' -^0 -^0 

xS{u;ok+i^ek' - t^p)aofc&efc,|0)o|0)e (1) 

for a co-linear pump. The sum extends over all possible 
wavevector modes. A polarizing beam splitter separates 
the ordinary and extraordinary beams into different spa- 
tial modes. The operators d and b represent these modes. 
Akk' is the phase mismatch, defined: kp — k — k' , tOp is 
the frequency of the pump laser, and the z integral ex- 
tends over the length of the crystal, fip and <jp are the 
central frequency and the spectral FWHM of the pump, 
respectively, e and o label the extraordinary and ordi- 
nary beams. The factor T) is defined as 41n(2). We are 
assuming the pump laser is Fourier-transform-limited. 

This state corresponds to squeezing just above thresh- 
old, such that mainly |0,0) and |1, 1) are produced. The 
|0,0) term is then dropped because it can not effect the 
two photon absorption process. 

Since the two beams are distinguishable after the po- 
larizing beam splitter, one of the beams must be subject 
to a polarization rotator in order for it to be made indis- 
tinguishable from the other. The beams must be indistin- 
guishable so that the |1, f) states interfere destructively 
and produce |2 :: 0) 

To study the two photon absorption probability of this 
state we utilize the well confirmed [20l[2l| Eqs. (2.15) 
and (2.16) from MoUow 



G(2) = {^iJ\E^-\t['^)E^-^t'2'^)E<-+\ti'^)EM{t2'^)\^i).{4) 

The primed and unprimed time variables represent the 
possibility of the biphoton traveling via two different 
paths of different lengths, as is the case in an interfer- 
ometric setup, and d labels the times as being detection 
times. (t) is the positive-frequency electric field op- 
erator defined by 

Sj kj 

xe-'"^.*^'a,^,^.(0). (5) 

Where the approximation e^*'^ '" « 1 has been made and 
s denotes either horizontal or vertical polarization. The 
time i'' = is defined as the time the photon is created. 
Note that = E'--'> . Also, f]/, and af are the central 

frequency and the FWHM of the filter in a specific arm, 
respectively. 

Now, by inserting a complete set of number states in 
the correlation function, and observing that all but the 
|0)(0| term will cancel, we can rewrite Eq. ([4]) as 

G(2)= (^^J\E^-\t['')E^-\t'2'')\0) 
x(0|£:(+)(ti'^)£^(+)(t2'')|1') 
~A{t'/,t'2'')*A{h'',t2'') 
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The above equation defines the biphoton ampHtude 
A{ti''-,t2'^). The expressions for A for SPDC were first 
calculated by Keller and Rubin in Rcf. 2^ and elabo- 
rated upon in Refs . [25l. l26l. l27t . We follow their calcula- 
tions somewhat closely. For more details see Appendix 
A. The resuh is A{ti,tf} = A{ti,4) + A{4,ti), where 
the script A represents the biphoton amplitude for down 
conversion followed by a polarizing beam splitter. The 
italic A is the biphoton amplitude for our setup. The 
result of our calculation for A is 



-i(Oet2''+noti'') 

UeUoUpaeCro<7p 
^ 7= 

X [Erf (T) - Erf(r + /)] . 



(6) 



Where 

r = 

I = 

And 



LU 



2UeUoUp^Via,^ + ao^ + ap^ 



Pu = Up{Ue - Uo)cTe<Jo, Eu = UeiUp - Uo)crp(To 
Ou = Uo{Ue - Up)<Jeap, = Pu^ + Eu^ + Ou^. 

The error function is commonly defined as 



Erf(.T) 



dye ^ 



The terms Ue, Uo, and Up represent the group velocities 
of the extraordinary, ordinary, and pump beams, respec- 
tively. The factors ae and do are the bandwidths of the 
filters in the arms of the interferometer. The term ap 
is the bandwidth of the pump. The factors fie and Qq 
arc the central frequencies of the extraordinary and or- 
dinary beams. For now the normalization constant has 
been left off for the sake of simplicity as it will only effect 
the height of the amplitude, not its overall shape. 

Hence we give the general biphoton amplitude in its 
most general form, which is a new result. 

In order to evaluate this expression we must find the 
group velocities of the pump, ordinary and extraordinary 
beams inside of the x^^'' crystal. The index of refraction 
as a function of wavelength can be found using the Sell- 
meier equations. Below are the Sellmeier equations for 
/3-Barium Borate (BBO) [H] 



no{Xo) - W2.7359- 



0.01878 
A2 - 0.01822 



0.0135A2 



ne(Ae) = J2.3753 



0.01224 
A2 - 0.01667 



0.01516A2. 



Where the wavelength is given in micro-meters. In type- 
II down conversion in BBO the pump beam experiences 
the same index of refraction as the extraordinary beam. 
It is also important to consider the angle the beams form 
with respect to the optic axis of the crystal. Since we are 
investigating the degenerate co-linear case with planar 
phase matching the optic axis must be set to be 42.4° off 
the pump beam's direction of propagation for a 400nm 
pump (the ordinary and extraordinary beams are co- 
linear with the pump and selected with a pinhole down- 
stream) [26j . For the pump and extraordinary beams we 
must use the effective index of refraction, given by 



COS^((/)) 



sm 



where is the angle between the beam and the optic axis 
of the crystal. We can now calculate the group velocity 



,eff 



^e,o,p(Ae,o,p) 



:(Ae,o,p) 



Ae,o,p '5'^ 



ff 

e,o,e 



(Ae,o,p) 



And tIq* ~ Uo- We then find for the degenerate case 
i^o = ^e) for Ap = 400nmthat Uo{^o) ^ 1.781 x lO^m/s, 
Uplilp) = 1.756 X lO^m/s and Ueifle) = 1.845 x lO^m/s. 

Fig. [3] contains a contour plot of the absolute value 
of the biphoton amplitude for our setup (A). This graph 
displays an interesting splitting which is symmetric about 
a line given by if — i2 = 0. Each point on this line rep- 
resents a different average arrival time of the biphoton 
{tf + t2)/2, a line drawn perpendicular to this line of 
symmetry represents another axis. That axis defines an 
entanglement time (the temporal distance between the 
two photons) as tf — f^- The symmetric splitting repre- 
sents the fact that two photons generated far away from 
the exit-surface of the crystal will drift apart in time. 
Thus, as the average arrival time increases (a delay being 
indicative of more time spent in the crystal) the photons 
drift apart. The symmetry is a result of the interferom- 
eter scrambling the information corresponding to which 
photon took which path. So for a set average arrival time 
there is an equal probability that the e photon will arrive 
first or that the o one will. 

We can check Eq. ^ by taking limits and comparing 
to known formulae. Taking the limit of the biphoton 
amplitude as <Jo and Ue go to infinity (the case of no 
filtering) and keeping in mind that one of the definitions 
of the Heaviside step function is. 



H{x) = lim -(1 + Erf(fcj;)). 

k — ^oo 2 



(7) 



We are left with 
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FIG. 3: Contour plot of the absolute value of the biphoton 
amplitude for our setup |^|. ae ~ cfo — (Jp = and 
L = 1.5cm. tf and t2 are in units of 10^^^ sec 



where 



A 
B 

J 



2UeUoVV 
L{Ue - Uo) 



2UeUoV-D 

[tfUoiUe - Up) + 4UeiUp ~ Uo)\ " 



V{Ue - UofUl 



Utilizing the identity 



Rcct \ -\ = U [x + - \ - H [x - - 



where 



Rect(a;) 



if > 0.5, 
\ if |a;|=0.5, 

1 if \x\ < 0.5. 



we obtain the expression 



yl(ti^,t2 ) = -2e-'"p-^Rect 



Which is equivalent to the expression given by Kim, et 
al., in Rcf. [11]. 



Eq. ([2]) given together with Eq. (fTT)) gives all the 
necessary information for calculating the two-photon ab- 
sorption probability for NOON states of = 2. Details 
of this calculation are given in Appendix B. The result is 



P2 = 



dvf 



P J-cx> 



|Erf(£:z./)-Erf(£z./-£)r 



1+4 



(8) 



Where 



1 



U 



V 



U2 

£ 
C 

C" 



1 

ul: 
1 1 

ITo'ul' 

iVV{U2(tI + Ueial + al)) 



LUi(ye^o 



2^V{al+al)' 



— V^)V27r7l? 

U2 



(9) 



And is the FWHM of the final state. The atomic re- 
sponse function g{uj) can be moved outside the integral 
for the light sources we consider (sec Appendix B for 
details). Unfortunately this integral is intractable an- 
alytically, we therefore perform the integration over Vf 
numerically. 



B. Numerical Calculation of Type-II Absorption 
Rate 



We wish to use Eqs. (|2l8l9p to calculate the relative 
absorption rate for type-II down conversion. After this is 
done, we can use the information to maximize absorption 
by adjusting the available spectral parameters. 



1. Pulse-Pumped Case 

We start with the case of a pulsed pump. Since the 
two-photon absorption probability is a complicated quan- 
tity with five adjustable parameters {ae^ ao, ap, L, and 
K/), we present several graphs to help elucidate the struc- 
ture of this mathematical object. All of the graphs 
are scaled so that the maximal absorption probability 
- within the parameters we consider - is set to one. 
Thus P2 = I does not represent an absorption proba- 
bility of unity. The normalization is consistent across all 
the graphs so that they may be compared to one another 
on the same scale. 
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FIG. 4: A plot of the scaled two-photon absorption probabil- 
ity for the realistic |2 :: 0) state as a function of the logs (base 
10) of the bandwidths of the filters in the arms of the interfer- 
ometer (in Hz). L = 2.3mm, Op = lO^^Hz, and k/ = 10^'* Hz. 



Fig. |3]is a plot of the two-photon absorption as a func- 
tion of the base-ten logarithms of the filter bandwidths. 
We sec that a wider filter is preferable. This is because if 
the spectrum of the light is wide then the temporal dis- 
tribution will be narrow, increasing the probability that 
the photons will arrive close together, triggering a two- 
photon absorption event. 

Fig. [5] is a plot of the two-photon absorption probabil- 
ity as a function of the logarithm base ten of the band- 
width of the pump for several different atomic linewidths. 
Clearly a broader atomic linewidth will result in better 
absorption, as this will allow a greater range of frequency 
pairs to be absorbed. For the pump beam a smaller band- 
width is preferable. Once the bandwidth of the pump 
approaches that of the filters, or the linewidth of the 
atom, the absorption drop off is dramatic. Furthermore 
a narrower pump bandwidth constrains the frequencies 
of the daughter photons, increasing the probability that 
the sum of their energies will be resonant. 

Fig. [HI is a plot of the two-photon absorption probabil- 
ity of a 1 2 :: 0) pulse as a function of crystal length for 
two separate cases. This is perhaps the most interesting 
of the representations of Eq. ([5]). The plot indicates the 
feature that for each setting of tio, fe, (Jp^ and Kf there is 
an optimal crystal length that maximizes the probability 
of two-photon absorption. 

As light travels through the crystal the dispersive na- 
ture of the medium causes the two photons to drift apart 
in time. This effect results in a drop off in absorption as 
length increases. Conversely, if the crystal length is very 
short, the spectrum of the light will be very broad, and 
the filters will strongly limit the amount of light that can 
reach the absorber. For broad filters, any deviation from 
the optimal length causes a dramatic drop off in absorp- 
tion. Information of this type would be useful when we 
design a two-photon quantum lithography experiment. A 



FIG. 5: A plot of the scaled two photon absorption probability 
for the realistic |2 :: 0) state as a function of the log (base 
10) of the bandwidth of the pump (in Hz) for three separate 
settings of the width of the final state of the absorber (in Hz) . 
L — 2.3mm and Oe = (Jo = lO^^Hz. 
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FIG. 6: A plot of the scaled two photon absorption probability 
for the realistic |2 :: 0) state as a function of the length of 
the crystal for two different settings of the filters (in Hz). 
Op = lO^^Hz and K/ = lO^Hz 



crystal cut to a specific length, for a given setup, would 
have the potential to enhance the two-photon absorption 
properties of the generated light. 

For example, take the case we have defined to be P2 = 
1: Here we have chosen, (Je = (Jo = 10^'^Hz, dp = lO^Hz, 
L = 2.3 mm, and k/ = lO^^'Hz. For the same absorber 
(k/ the same), ae ~ (Jo = lO^^Hz, ap = 10^'^Hz, and L = 
2 cm, the absorption probability is is P2 = 2.07 x 10~^. 

We may obtain completely analytical results for the 
case where there are no filters in the arms of the inter- 
ferometer. Using Eq. ([7]), Eq. ([8]) becomes 



lim 



2L(Jr, 



-Ko 



(10) 



so long as i ^ 0. The function Kn is the modified Bessel 



Where 



2.xl0'2 4.xl0'2 6.xl0'2 8.xl0'2 1.x 10" 
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FIG. 7: A plot of the scaled two-photon absorption proba- 
bility for the 1 2 :: 0) state (in the limit of no filtering) as a 
function of the bandwidth of the pump. We indicate three 
separate settings of the width of the absorber. We also set 
L = 1 mm. The scaling is the same as on the other pulse 
pumped plots. 



function of the second kind. Features of this function 
include a simple dependence on L and a modified Cp 
dependence, as shown in Fig. [71 Interestingly, for ab- 
sorbers with very wide final states, the two photon ab- 
sorption probability becomes almost independent of the 
bandwidth of the pump. 

In all cases, given a particular Kf, a balance must be 
struck between the related quantities of temporal and 
spectral correlation. Tight temporal correlation increases 
the probability of absorption. However this necessitates 
broad spectral distributions which have a lower probabil- 
ity of matching the correct energy for transition. 



2. Comparison to Coherent Light 

We can now tell how varying the spectral properties 
of 1 2 :: 0) states will effect the absorption probability. 
However we have discovered nothing about how |2 :: 0) 
light's absorption properties compare to other, more fa- 
miliar, states of light. We would like to say wether |2 :: 0) 
has any benefits or disadvantages when compared to, for 
example, coherent light. Coherent light has the advan- 
tage of being relatively well understood (or at least very 
extensively studied). Thus we shall attempt to find a 
state of coherent light which will serve as a fair com- 
parison to 1 2 :: 0), as produced in the fashion presented 
in this paper. We would like this fair state to have the 
same spectral profile as |2 :: 0), but to lack the unique 
properties that the NOON state possesses: momentum 
entanglement and a high degree of temporal correlation. 
For details on the derivation of this state see Appendix 
C. The result for the two photon absorption probability 
is 



_ Lupap 
" ~ V2P ■ 

(12) 

and We = jj -jj- and Uo = -jj jj-- Here, / is the 

intensity of the light in each arm of the interferometer 
(i.e. 21 would be the total amount of light in the device). 

This fair-comparison light can be though of as a very 
spectrally broad coherent source incident on a magic "fil- 
ter" , which imposes on it the same spectral profile as the 
1 2 :: 0) state. Two spatial modes of this light are then 
used as the input into an interferometer without filters 
in the arms. The two-photon absorber is placed at the 
far end. 

The error functions in the the above coherent-state ex- 
pression diverge much more rapidly than for the |2 :: 0) 
case. Thus it is only possible to compute for rela- 
tively narrow arguments of the integral over Vf, so we 
must consider a non-maximal setup as a test case. So for 
CTp = lO^Hz, Kf = 10^"Hz, and L = 1 cm we have 

^ 5.65 X 10"^P2/^ 

where has been calculated using Eq. ([7]). 

If the intesities of the two kinds of light are set to 
be equal (/ = 1 in the above equation) then |2 :: 0) is 
absorbed with a much higher probability than coherent 
light. This is directly a result of the frequency entan- 
glement. Maximal absorption probability occurs when 
the photon energies add, such that the total is the same 
as the center of the final level. For SPDC it is insured 
that this condition will be closely met (exactly met in 
the case of a cw pump and a very thin crystal). However 
for the coherent analogue there is no such requirement, in 
fact the probability of two randomly chosen photons from 
the coherent pulses adding up to the resonance energy is 
miniscule. Even so, the above equation may be overly op- 
timistic about the absorption properties of coherent light. 
Recall that we were forced to consider narrow-band fields 
to make the calculation numerically tractable. The ef- 
fect that causes coherent light to have a poor two-photon 
absorption probability will be worsened when the band- 
width is broader; a regime in which the |2 :: 0) states' 
absorption probability improves. 
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These positive results for the muhiphoton absorption 
of entangled light is in seeming contradiction to the 
results of section II, which shows that coherent states 
should fair better than \N :: 0) states, where TV-photon 
absorption is concerned. However that treatment does 
not take into account the spectral properties of highly 
quantum- mechanical states of light, which arc likely to 
have high degrees of temporal correlation. 

In a paper by Tsang [lO| , it was shown that the spatial 
properties of NOON states cause them to have poor N- 
photon absorption rates. The relationship found in that 
paper is that NOON states have absorption rates that are 
lower than a classical analogue (in his case a monochro- 
matic Fock state incident on a beam splitter) by a factor 
of l/2^~^, due to spatial considerations independent of 
temporal correlations. He leaves open the question of 
whether time-domain effects can compensate for this ef- 
fect. Given our above results it seems likely that the 
answer to this question is yes. 

However, the advantage entangled light enjoys is offset 
by the fact that SPDC has an extremely small intensity. 
Even relatively high production efficiencies are only ap- 
proximately one pair per 10^^ incident photons (sec for 
example Ref. [1^). Furthermore it has been shown that 
correlated two photon absorption (absorption from two 
photons of the same pair as opposed to absorption from 
two photons of different pairs) from entangled light dom- 
inates only when the intensity is small [30|. This makes 
it even more important to improve absorption rates. 



where 

In the second equality we have taken the continuous limit, 
changed variables from momentum to frequency, per- 
formed one of these integrals to eliminate the delta func- 
tion, and then changed variables again to {i^el ~ \i^o\ = 
We are able to do this due to conservation of energy for 
a cw pump: once the momentum of one photon is chosen 
the momentum of the other one is determined. Also the 
phase mismatch becomes 

We now write in analogy to Eq. pT]) 

g-i(a„ii''+aet2'') 
A{h',t/U = jjjj-^ 

X / dzT{v ^ [is'' - ti''])Y{v) 
Jo 

The solution for A is easily found 



3. Continuous-Wave-Pumped Case 

Let us now investigate the case of a continuous-wave 
pump. The expressions derived in the previous section 
only make sense when considering a biphoton pulse. They 
give the absorption probability of a biphoton of finite ex- 
tent being absorbed as it passes an atom or other ab- 
sorber. Note that since the pulses are Gaussian-like they 
are not strictly finite, but they are so closely temporally 
correlated that we may treat them as such. However the 
continuous wave case is stationary. We must now speak 
of a two photon absorption rate. 

First we recalculate the biphoton amplitude. Take Eq. 
p4p . for a cw pump we set cjp = and ujk^ + i^k^ — ^p- 
We can then rewrite this equation as 



^(ii^t2'')cw= ^e^"'(■ 



d ± d\ 



-i(0„ti''+0^t2'') 



Aiti^M 



Ue-Uo 



Erf 



{h'-t^'^y 



2VV 



Erf cr 



{ti"^ - t2'^) - Lu 
2VV 



The absorption rate for stationary states is given by 
Eq. (3.17a) in Ref. [H 



W2 = 2 



.9 ( 



dte^^^pt-'^f\t\G^^\-t,-t;t,t), 



where k/ is the width of the final state of the absorber. 
We have for the second-order correlation function, 



G^'^\-t,~t-t,t) 



2e 



[Ue - Uof 



Erf 



Lua 
2VV 



and thus 



Jo 

uu / / 

UeUo J-oc Jo 



xe 



W2 



2|g(i»p)|' 



Erf 



Lua 
2VV 



The modulus squared of the error function is the only 
part that is dependent on the spectral properties of the 
field, so we fold the overall constants into the atomic 
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FIG. 8: Plot of the two-photon absorption rate from a cw 
pumped crystal, as a function of crystal length, for four sep- 
arate settings of the filters (in Hz). 




FIG. 9: Plot of the two-photon absorption rate from a cw- 
pumped crystal, as a function of the bandwidths of the filters 
in the ordinary and extraordinary arms. Here L — 6 mm. 

response function, which we wiU then ignore. We then 
obtain the relatively simple result 



So for a cw pump we have generally the same results 
as in the pulse-pumped case, with absorption improving 
as bandwidth increases. We also again observe that the 
graph of W2 as a function of L displays peak values. 

IV. CONCLUSION 

First, we analyzed the multiphoton absorption proba- 
bility of states of the form :: 0) in the ideal case (with 
no spectral information). We found that the absorption 
probability scales poorly with N when compared to co- 
herent states, but well when compared to Fock states. 

Second, we considered the case of a realistic |2 :: 0) 
state produced by type-11 spontaneous parametric down 
conversion, a polarization rotator and a beam splitter. 
We found that generally the two photon absorption prop- 
erties of this states arc highly dependent on the specific 
optical setup used. Using numerical methods it is pos- 
sible to analyze the absorption probability as a function 
of the realistic parameters of the experiment. We can 
then adjust the optical setup of our model, including the 
length of the crystal, the bandwidths of the filters in the 
extraordinary and ordinary beams, and the pulse length 
of the pump. Running a maximization procedure over 
these variables it is possible to find the setup, which op- 
timizes two-photon absorption. The difference between 
the optimal setup and a slight deviation may be dramatic. 

This research constitutes a proof of principle of the idea 
that the poor production rates and detrimental spatial 
effects of highly quantum mechanical states of light may 
be mitigated by improving the absorption rates through 
spectral means. Though we consider in detail only the 
1 2 :: 0) case, it is likely that once methods of developing 
NOON states of higher N are developed, similar methods 
may be applied to improve their N photon absorption 
properties as well. 

Also, in a broader sense, similar techniques may be 
applied to increase the multiphoton absorption properties 
of any desirable state of light. All that is required is the 
quantum mechanical state vector. 



where the factor of 1/L comes from the state normaliza- 
tion. Unlike for the general case, the cw expression is 
fully analytical. Figs. [8] and [9] are plots of the absorption 
rate as a function of crystal length and filter bandwidths. 
The graphs have been normalized such that the greatest 
absorption rate for the range of parameters we consider 
is set to W2 — 1. 
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Appendix A: Calculation of the Biphoton Amplitude 

We start with a field at the detector, given by Eq. 
Using the standard 50:50 beam splitter transformation 
we write out the electric field operator at the crystal 



^/2 



Where constants have been subsumed into the overall 
factor of ^. The polynomials ujkj have been ignored be- 
cause they vary slowly when compared to the exponential 
terms. The operators a and b are the annihilation opera- 
tors acting on the first and second modes. The two mode 
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annihilation operators are used to signify that the amph- 
tude is dependent on the fields in both spatial modes. 
T = — l/c where I is the distance between the crystal 
and the detector. We only consider the case where the 
interferometer is path balanced. This implies we calcu- 
late the absorption properties of the central fringe. We 
assume that the absorption properties of the other fringes 
will behave similarly. 

Now using the above and Eqs. (HI 21) we can find an 
expression for the biphoton amplitude of our setup 



(14) 



where A is the biphoton amplitude for just the output of 
BBO and filters 



xe 

/•^ - 
X / dze 
Jo 



xe ' e 

Where Ue = y- and Uo = jj y-- Note that 

dk ~ — jjduj = jjdh', so that up to a constant we 

may switch integration variables between momentum and 
frequency. (Actually the t/'s are frequency dependent, 
however they do not vary significantly over the band- 
width of the field.) Here we diverge from Kim and Shih 
in that we integrate over Vo and i/g instead of = j/q — i^e, 
and i^p and we do not take the filter bandwidths to be 
equivalent. Now 



Aih^^h^y- 



UeUo J J JO 



We have dropped the factor of e^*'/"^ as it just introduces 
an overall phase. The central frequencies of the filters 
have been chosen to be the same as for the e and o rays. 
Now, define the following variables 



(15) 



Due to the delta function in Eq. flo + fie ^ and 
i^o + i^e = Vp. Now, Taylor scries expand the wavevector 
out to second order 



kivj) = k{nj) + vj 

= k{VLj) 



dk{vj 



dujj 



Where j = e,o,p and U is the group velocity and U 
doj/dk. Taking into account that kp{Q.p) = k^i^e) 
ko{flo) The phase mismatch can now be rewritten as 



UpiQp) Uoino) C/e(r!e) 

Vq + _ Vq Ve 

UpiVlp) ^ Uo{no) ^ Ueifle)' 



(16) 



Taking the continuous limit of Eq. (|14p and utilizing Eqs. 
rSl) and ITBl) wc obtain 



We have the unitary Fourier transform over i/q, Ve (up 
to a factor of V2tt), and the integral over z given by 



Jo 



T(z/o,j/e,z)= e 



xe 



iz{Uel'e--\-Uoyo) 



(17) 



Where we have defined the biphoton kernel T(i'o, VeTz). 
In the interest of clarity the Fourier transforms have been 
written as operators 



T{x y) 



J- 



dxe-'^'y 



Using Mathematica"^^ , we perform these operations in 
the above order. Note that another ordering will lead to 
the problem becoming intractable. The result is given by 
Eq.®. For the sake of simplicity in most of our calcula- 
tions we will take ^ = 1. 



State Normalization 



The state needs to be normalized. Take again 
Eq. ([T]) where the integral over ujp has been performed 
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Now wc take the inner product 



alal\Q), 



(*l*> 



4C^ 



O J — OO ^ — CX3 



Sin^ [z/gWe + i^o^o]) 
X 2 ■ 

Where wc have taken the continuous hmit and switched 
integration variables. It is required that = 1. We 

make another change of variables: p = i^^ + i^o and q = 
Vf.Uf. + VqUo, as long as ^ Uo this is a well defined 



one-to-one transformation with a Jacobian J 
This allows us to separate the integral. Thus 



4C2 



C/eC/o(f/e - Uo) 



dp e K'pJ dq- ' 



2C2 



■n '^ La. 



p 



UeUo{Ue-Uo) V2V 



So 



This factor is added onto Eq. ([8]). Using a similar pro- 
cedure we obtain for the continuous wave case 



UeUoiU, ~ Uo) 



2tiL 



which we append to Eq. 



Appendix B: Calculation of the two photon 
absorption probability for NOON states 

Noting that Eq. ^ represents a unitary transform 
back from time space into frequency space we may write 
the Z function as 



Z{uJo,UJe) = 2T:T{h'' UJoTHh'' UJ,)*A{ti\t2''). 

Therefore 



Jo 

= C / dz:F{t2'' V,)*P{U, ^ t2^)^(il'' ^ i^o)*^(i^o ^ tl'')T(l^o, J/e, z) 

Jo 

+ C' I dzJ^ih'^ ^ VoTT[v, ^ tY')T{t2'' ^ VeYTiVo ^ i2'')T(l/o, i^e, ^) 

Jo 

= C' f dzT{yo,Ve.z) + C' I dzT{u,,iyo,z). (18) 
Jo Jo 

I 

In the second equality we have made a substitution using T(t/e, 1^0, z). And 
Eq. HH). In the third equality Eq. p7|) was used. In the 
fourth equality the overall phases have combined with the 



Fourier transform operators to change their output vari- ^ fUi 

ables. So the Fourier transforms from the biphoton am- C = 27ri(87rP)* 4 / -j— 

plitude exactly cancel with the Fourier transforms from V ^ 

the definition of the spectral correlation function, leav- 
ing a relatively simple result. Note that T{uo,ye,z) + '^'^^^'''^ ^2 is defined in the main text. Returning to 

Eq. ^ we make the assumption that g{io) « g(}^Sly). 
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To see why, we need to investigate the definition of the 
atomic response function, Eq. (3.8) from MoUow 

5(w) = tJ-'^PSiPjo ^— 1 — • 

Where ^ is a constant representative of the absorber, j la- 
bels the possible intermediate (virtual) levels, pj^ and pyj 
are the momentum matrix elements of the electron mak- 
ing a transition between the initial and jth states, and 
the jth and final states respectively. Kj is the linewidth 
of the intermediate level. Each term in the series repre- 
sents the two photon absorption process proceeding via 
a different intermediate level transition. We assume that 
the central frequency of the field is one half the resonant 
frequency of the final state. Note that each term of g 
is highly peaked around the characteristic frequency of 
the level, which will be far detuned from the central fre- 
quency of the field. Unless the bandwidth of the field is 
on the order of the transition frequency, the value of the 
atomic response function will not change much as uj is 
varied over the bandwidth of the incident light. There- 
fore our assumption is justified. Note that this is the only 
assumption we make about bandwidths in this paper. It 
is worth noting that there has been much work done to- 
wards enginee ring materials which have large two photon 
cross sections [Ull^], so as light sources are improving, 
so are absorbers. 

So, using Eq. and Eq. ^ we have 



the pump beam being absorbed. Furthermore it implies 
that this frequency will be on resonant with the final 
level. Realistically, all frequencies will have a chance to 
be absorbed, and only Vf = {) will be resonant. To correct 
for this we can average the function over a lorentzian line 
shape with its peak at = and a height of 1. The full 
equation is given in the main text by Eq. ([8|). 



Appendix C: The Fair Comparison Coherent State 

Again, our |2 :: 0) state is produced from a photon pair 
from a type II spontaneous parametric down conversion 
event, incident on two filters and a symmetric beam split- 
ter. We can examine |^), Eq. ([1]), and see that - given 
the assumptions in Appendix A - this state has a spectral 
profile (the relative probability amplitude of the photons 
having frequencies z/g and Vq) given by 

Jo 

1 4") represents two entangled down conversion modes. To 
remove the non-classical nature of this light we project 
out one of the modes. We arbitrarily choose the ordinary 
mode. The remaining state will have the spectral profile 
we desire for one mode of our fair coherent state (Fa), 
that is, 



A = 



C'g (^np) f dv [ dz 
X [T{j,f - V, V, z) + T[y, lyf - z)] \\ (19) 



The importance of the atomic response function being 
separated from the integral should be noted. Essentially 
this means that the spectral properties of the light may 
be considered apart from the structure of the absorber. 
The remaining integral is not difficult to perform 



dz 



dvT{vf — z) 



X [Erf {£i^p) - Erf {Siy,, - £)] . 

Where t/2, We, £, and C are defined in the main text in 
Eq. The integral over T{i', Vf — v, z) differs only by 
a constant factor. 

However, the utility of this expression is limited. It 
only describes the probability of a single frequency from 



dVo{luA^) = \/-^^Fa,,(i/e)aL'|0>e, 

-°° k' 

"'0 



where = ^/Up — 1/Ue and Uo ~ ^/Up — l/Uo- This 
will represent one mode of the "fair comparison" coher- 
ent light. The normalization will be discussed shortly. 
The other mode will have a separate, but identical, spec- 
tral profile. Both modes will be mixed with a symmetric 
beam splitter. 

In order to assimilate this spectral information into 
a coherent state we identify the spectral profile with a 
thus: Uk = Fa^{i'k). 

The beam splitter transformation operating on the dis- 
placement operators which create the two mode coherent 
state, |a)i|/3)2 produces an output of |^i±|^)3|^)4 . 
Where 1 and 2 label the input modes and 3 and 4 label 
the output modes. 

So now we can write the temporal correlation function 
for the fair coherent state 
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^ ^^^^^ 
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4f ,y„ 1/^ 



a + if3 



V2 / 3^ 



^ I it, 



Where the factors 7 are the normahzations in the indexed 
modes. The central frequencies (f2) have been taken to 
be the same. The filters in the arms ~ Uo have been 
removed in order to make the calculation which follows 
mathematically tractable. The frequencies have been la- 
beled a,6,c, and d. As before ^ represents the constants 
associated with the electric field operators. We shall 
eventually take the two two photon absorption proba- 
bilities in ratio, the ^'s will cancel exactly. In light of 
this we simply set ^ = 1. We need to normalize such 
that the information about the intensity of the coherent 
light in each arm is contained in the 7's. 



spectral correlation function 



5(2) = 2^^\^,\^\^,\^[F:[v'^)-^F*p{u\)\ 
X [F*p{u'2)-iF:{u'2)]. 

Now we can write the two-photon absorption probability 
Eq. ^ as 
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= I7P / dv' / dzi / dz2 

J-co Jo Jo 

= ["^dz, r dz,Siz,-Z2)e-^^^i+^'^^ 

Uo ~ We Jo Jo 

The double integral is tractable yielding the result 



Using the same procedure as in Eq. psp we obtain the 
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-00 ""o — Ue Jq Jq 
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x|Erf(£„z//) - Erf(£„j// + £„) 



Where £a and Ca are defined in the main text. In the 
interest of simplicity the initial coherent states have been 
taken to be equivalent. As in Appendix B we must 
now take the integral of the above expression, times a 
lorcntzian line shape, over Vf. This new function, P^f , is 
given in the main text by Eq. pip . 



